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Abstract. These notes develop aspects of perturbation theory of matrices related 
to so-called diagonalisation schemes. Primary focus is on constructive tools to derive 
asymptotic expansions for small/large parameters of eigenvalues and eigenprojections 
of families of matrices depending upon real/complex parameters. 

__ ■ Applications of the schemes in different frameworks are also discussed and references 

j.-^ ' to further applications given. 
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H . 1. Introduction 

In various situations one is confronted with the question of determining eigenvalues 
and eigenprojections of matrices uniform in certain parameters. The available tools for 
answering this question range from hand calculation (which is fine if matrices are small, 
i.e., have size 2x2) to strong abstract methods of perturbation theory providing analytic 
formulas (but no concrete answer in most cases). In this note we want to collect results 
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on a particularly useful constructive scheme for calculating asymptotic expansions for 
diagonalisers and therefore also eigenvalues and eigenprojections. The method is flexible 
enough to generalise to applications to systems of differential equations or matrices of 
pseudo-differential operators. 

The approaches presented here generalise those used by Taylor [TB], Yagdjian [23], Wang 
[T7], [18] and others, including the authors [12], [Ij. 

The paper is organised as follows. First we will recall in Section [2. II the well-known case 
of matrices with distinct eigenvalues and their small perturbations. In Section 12.21 we 
will present the general multi-step diagonalisation scheme, and its relations to uniformly 
diagonable matrices will be discussed in Section [231 Section [3] is devoted to applications 
of the schemes in different situations. 

2. Diagonalisation schemes 

2.1. Non-degenerate matrix families and the standard scheme. We start our 
presentation by recollecting some well-known aspects from perturbation theory in com- 
bination with a merely classical proof using a diagonalisation technique. Let for this 
A : M — > C™"^™ be a continuous matrix-valued function depending upon a real (or for 
some applications also complex) parameter p and assume that as /o — >• the family of 
matrices has a full asymptotic expansion 

(2.1) A{p)^ Ao + pAx + p^A2 + --- , p^O, 
i.e., 

(2.2) A{p) = Aq + pAi+p^A2 + --- + p^An + 0{p^+^), ^N, 

for certain (uniquely determined) matrices Ai G ([^rnxm^ ^^ ^j^g ^^^^ q£ complex param- 
eters this is just slightly weaker than assuming analyticity of A near p = 0; in the case 
of real parameters the assumption follows from smoothness oi A in p = 0. 

Now one might ask how eigenvalues and eigenprojections of the matrices A{p) depend 
on p as p ^ 0. It is well-known that continuity of A{p) implies that eigenvalues depend 
continuously on p, but in order to conclude more we have to make assumptions on A{p). 
It is reasonable to define a non-degenerate matrix family as one with distinct eigenvalues. 

Definition 2.1. We call A{p) non- degenerate in p = if Aq has m distinct eigenvalues. 

Theorem 2.1. Assume that A{p) is non- degenerate in p = 0. Then there exist uniformly 
bounded families of invertible matrices M[p) with uniformly bounded inverse having full 
asymptotic expansions as p — > and satisfying 

(2.3) A{p)M{p) - M(p)A(p) = O(p^), ViV, 
for a diagonal matrix A(p) . 

The diagonal matrix A(p) coincides to arbitrary order with the (continuous) diagonal 
matrix containing the eigenvalues of A{p). It has an asymptotic expansion and the 
corresponding coefficients are determined uniquely. Although the above statement looks 
rather weak, it is of particular interest for us. We will show how to construct the 
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diagonaliser M{p) using a recursion scheme and how the aforementioned uniform bounds 
and the asymptotic expansion of A(p) arise naturahy within the construction. 

The construction forms the core of the more involved multi-step scheme introduced later 
in Section 12.21 and is the key idea in the diagonalisation-based approaches of Yagdjian, 
Reissig and co-authors, see e.g., [23], [13], [9], [8]. For more detailed discussions on 
applications we refer to Section [3l 

Proof, (of Theorem I2.ip It is enough to prove that for any number N there exists an 
interval Im = {—^n^^n)-, a uniformly bounded and invertible matrix function Mi^{p), 
p € In, and a diagonal matrix Ajv(/)), p £ I^, such that 

(2.4) A{p)Mn{p) - Mn{p)An{p) = 0{p''+^) 

and M7v+i(p) — Mn{p) = 0(/7^+^). Then any two functions M{p), A{p) coinciding with 
all the Mj\[{p), A]\f{p) up to the corresponding orders satisfy the above theorem. But 
this just means we need to construct functions subject to asymptotic expansions, which 
is a standard argument of asymptotic analysis. 

We are going to construct M^ip) and Ajv(p) recursively, and so assume that they can 
be written as 

(2.5a) Mn{p) = Mo{I + pM^^^ + p^M^^) + • • • + p^M^^^), 

(2.5b) A^ (p) = Ao + pA(i) + • • • + p^ A W 

with certain coefficients Mq, Aq, M^^^ and A^-^^. It is clear by (f2lil) with N = that 
AqMq = MqAq and so Mq must be the diagonaliser of the main part Aq and Aq the 
corresponding diagonal matrix consisting of all eigenvalues of Aq. 

Step 1. We determine a diagonaliser Mq of Aq and set A{p) = Mq A{p)Mq. The new 
matrix family A{p) has the asymptotic expansion A{p) ~ AQ + pAi+p'^A2 + • • • as /> — > 
with coefficients Af^ = Mq Aj^Mq. 

Step 2. Assume we have already determined Mk-i{p) and A^-iip) for a certain number 
A: = 1, 2, . . .. Then we denote 

(2.6) B^''\p) = A{p)Mk-i{p) - iV4_i(p)Afc_i(p) = 0{p^). 

It is evident that this matrix family has a full asymptotic expansion as /) ^ 0, and it 
makes sense to set B^ ' = livUp-^Q p^ Mq B^ \p). We now construct the next coefficient 
matrices. For this we set A^'^) = diagi?^'^^ and take M^^> to be the solution of the 
Sylvester equation 

(2.7) [Ao, M(^)] + i?(^) - diag B^^^ = 

with vanishing diagonal entries. Equation (12. 7p is soluble for any B^ \ because Aq = 
diag(Ao,i, . . . , Ao,ra) has distinct entries, and its solution is given explicitly by 

f (sW) 

(2.8) (mW) = "V"^' i^J^ 

0, i = j. 
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That these are indeed the right choices for the matrices follows from 

5(^-+i)(p) = sW(p) + p^ (a{p)MoM^''^ - MoM^^^Akip) - A4_i(p)A(^)) 

(2.9) = p^Mo (^('^^ + [Ao, M^'^)] - A(^')) + 0(/+^) = 0{p''+^), 

by (|2.7|) and so the construction can proceed recursively. 

It remains to check invertibility of M]\f{p) on a suitable In- The matrix Mq is invertible; 
the second factor of the form I + 0{p) and tends to the identity as p ^ 0. Since the 
group of invertible matrices is open in C™^™, this implies invertibility of Mj\f{p) for 
sufficiently small /at- D 

Remark 2.1. Determining the diagonaliser Mq in Step 1 of the scheme is the only 'painful' 
step of the approach and requires considerable effort. Everything else is explicit and 
consists of operations acting on entries of matrices (in (12. 8p ) or matrix multiplications / 
additions (in (12.61)1. 



Remark 2.2. The intervals Im will in general shrink as A^ tends to infinity. This is 
different if A[p) is analytic in p near and the corresponding asymptotic series converge 
uniformly on a small interval (see e.g. [51 Chapter II]). In this case (j2.3|) implies that 
the error term is actually 0. 

For completeness we also mention the following spectral bound, which can be used to 
estimate the eigenvalues of A{p) in finitely many steps of the diagonalisation scheme. 



Proposition 2.2. Representation ()2.6p implies 

(2.10) dist( spec A(p), spec Afc_i(p)) < \\Mj:^^{p)B^''\p)\\ = 0(/). 

Proof. The diagonal matrix Afc_i(p) is normal, and so its resolvent satisfies the bound 

(2.11) ||(C-Afc_i(p))-i||< ^ 



dist(C,specAfc_i(p)) 

Thus for any C, £ C with dist(C,spec Afc_i(p)) > \\Mj^\{p)B^^\p)\\ the right-hand side 
of the resolvent identity 

(2.12) (C-Afc_i(/>))-'(C-M-_\(p)^(/,)M,„i(p)) 

= I-{C-K-i{p)y'M~\{p)B^'^\p) 

is invertible. But this implies invertibility of the left-hand side and so that C, belongs to 
the resolvent set of A[p). D 

Remark 2.3. The bound is almost optimal. If p is sufficiently small the right-hand 
side of the estimate involves essentially all entries of B^ \ while the next step of the 
diagonalisation brings in the diagonal entries of B^^' as new coefficients. So the estimate 
is sharp if B^^> has no exceptionally large off-diagonal entries. 

Another interesting consequence is that the constructed matrices M^ip) allow to ap- 
proximate the eigenprojections. To make this precise we consider the eigenvalue Xj{p) 
with its eigenprojection Px./p). It corresponds to the j-th diagonal entry of the matrix 
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Ai\f{p). For the diagonal matrix it is evident that the corresponding eigenprojection is 
just Bj (8> ej for the j-th basis vector Cj of the standard basis. The eigenprojection P\ [p) 
is determined by (Aj(p) — ^(/j))Px (p) = (and the fact that it is a projection, i.e., 

^Im =^a,(p))- 

Proposition 2.3. The projection M]^{p){ej®ej)M^ (p) approximates the eigenprojec- 
tion P\ (p), in the sense that 

(2.13) {\,{p) - A{p))MN{p){e, e,)M^\p) = 0{p''+') 
and 

(2.14) ||P,^,(,) - Mjv(p)(e, e,)M^i(p)|| = 0(p^+i). 



Proof. The first formula is a straightforward consequence of (j2.4p and Proposition! 
For the second one we recall that the eigenprojections of a matrix can be represented in 
terms of the eigenvalues by the product formula 

(2.15) Pa,(p) = l[(^^ip) - X,ip))-H\ip) - A{p)). 

Plugging in (j2.4p and (j2.10p for A{p) and Xj{p), respectively, the second statement 
follows. In both cases the constants in the estimates can be calculated explicitly. D 

2.2. Block-diagonalisation. The main objective of this note is to discuss how to gen- 
eralise the scheme from the proof of Theorem 12.11 to degenerate matrix functions and to 
replace the assumption of non-degeneracy by weaker assumptions, which are just enough 
to ensure the existence of an asymptotic expansion of eigenvalues and eigenprojections 
in entire powers of p. 

The main motivation for the approach goes back to Taylor [16j and Wang [T7], [18| 
and applications to the decoupling of hyperbolic-parabolic coupled systems. The precise 
construction is taken from the PhD thesis of the first author, [3]. 

A short comment on notation: for matrices appearing recursively in the scheme we use an 
upper index [j] to denote the level. Matrices which are final results of the consideration 
get lower indices according to their position and if these matrices are themselves families 
with asymptotic expansions we use upper indices (j) to denote their j-th term. 

2.2.1. Basic setting and diagonalisation modulo 0{p^). We assume that the matrix Aq 
is diagonable, and denote by Mq a diagonaliser of Aq, which arranges the eigenvalues in 
groups, i.e., 

(2.16) Mq-i^oMo = Ao = diag(Afl, • • • , A^) 

and \\ = \- \ i < j , implies A^ = A^ for all i < k < j. We apply Mq to the family 
A{p): 

(2.17) iM(p) = M^^A{p)Ma ~ Ao + pAf + p'Af + • • • , p ^ 0. 

We assume that Aq is not just a multiple of the identity and try to follow the standard 
scheme as far as possible. If we consider the Sylvester equation [Aq, X] = Y , we see that 
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this can only be solved if {Y)ij = whenever Aj = Xj. More formally, we introduce the 
notation Hq = (tti, . . . , Tr^), vri + • • • + vTfc = m, for the partition of eigenvalues, and write 

i ~no J if -^i ^j ■ Furthermore, b-diagnp selects the corresponding diagonal blocks 
out of a matrix / constructs a block-diagonal matrix of this form. 

Let Ai = b-diagng ^i ■ Then we can solve the Sylvester equation 

(2.18) [Ao,i^fU4°'-Ai = 0, 
and obtain as possible solution (cf. (j2.8p ) 

(4°') 

(2.19) «%,. 



0, i ~no J- 






Furthermore, 

(2.20) iM(p)(/ + pkS'^) - (/ + p4°1)(Ao + pAi) = /,(iM + [Aci^f^] - Ai) + ©(p^) 

is of order p^. In order to get the matrix diagonalised modulo 0{p^) it remains to 
diagonalise Ai. For this we observe the following: if Ai is diagonable, its diagonaliser 
lives in the subspaces corresponding to the blocks / the partition Ho. On all these 
subspaces the matrix Aq is a multiple of the identity (by definition of the partition). 
This implies that Aq is invariant under all diagonalisers of Ai, provided they exist. 

Thus we assume that Ai is diagonable, and denote by Mi a diagonaliser of Ai, arranging 
the eigenvalues into groups within the the partition Ho. Then 

(2.21) Mo(I + pi^i°^)Mi 

diagonalises A[p) modulo 0{p^) for small values of p. Hence we constructed 

iW(p) = Aff 1(1 + pi^f])-^iM(/,)(/ + /54°l)Mi 

(2.22) ~ Ao + Ml + P^A^2 + • • • ' P ^ 0. 

Assumptions we had to make were the diagonability of the two matrices Aq and Ai = 
b-diagngMo-^iMo. 

Remark 2.4. Instead of diagonalising Ai and stopping the procedure we could also apply 
the construction iteratively and obtain a block-diagonaliser 

(2.23) Mi{p) ~ / + pM(^) + p^M^^) + • • • , 

(with M^^' = i^Q and the further terms obtained by a similar procedure to the standard 
scheme) such that M^^{p)A{p)Mi{p) is Ilo-block-diagonal modulo f]j^O{p^). 

2.2.2. The iterative scheme. We assume we already applied k steps to diagonalise the 
given family modulo 0{p^^), i.e., we assume we are given 

(2.24) iW (p) ^ Ao + pAi + • • • + p'^Ak + Z+'ifli + • • • , P^O, 

and in particular that the eigenvalues of A^, are arranged into groups within the partition 
Hk-i- Associated to A^ we have a new partition 11^, which is a sub-partition of II^.i, 
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and a corresponding equivalence relation ~nj. • Denoting the entries of A^/ as A^ , . . . , 
Am this means 

(2.25) i^n.j ^^ yk' < k : Xf'^ = Xf'\ 

~\k] 

Our strategy now is as follows. The matrix Aj^'-^ is a full matrix, which does not obey 
the special block structure determined by the filtration 11^ of partitions, but we can 
construct step by step diagonaliser which eliminate off-diagonal terms modulo 0{p^~^'^) 
with respect to the 11^/, k' < k. This can be done in A; + 1 sub-steps. 

\k 01 ilk] \k] 

For the 0-th step let Af^^^ = b-diagug ^fc+i and Kq be the solution to the Sylvester 
equation 

(2.26) [Ao,Kn + 4'|i-4+'? = 0. 



Then by construction 



(2.27) iW(p)(/ + p'^+^<') - iI + p'^'K^)iT.f^^^+P'^'4w 

is of order 0{p^~^'^), the last term on the right of Ilo-block-diagonal structure. 

In the ^-th sub-step, i = 1, . . . ,k, we denote ^^_f_| = b-diag^^ ^L+i > ^^^ define K^ 
to be a n£_i-block-diagonal solution of 



(2.28) [A,,i^fi]+4Yr"-4Yl = o- 



Wi 



Again by construction and the commutation property [A£i,K^ ] = for ^ < ^ (because 

\k] 

Kg acts only on the invariant subspaces related to n£_i), it follows that 

(2.29) (X:p'A,+p^+i4Yr'')(^+/''"''~'^f)-(^+p'^'~'^f')(E/^'^^-+/''^'4+i) 

i=o j=o 

is of order 0{p^'^'^). The last term on the right has Il^-block-structure. 

Finally we obtain a block-diagonalisation up to Il/j-structure with last remaining term 

T \k k] 

Ak-\-i = ^fc_Li • If we assume that this matrix is diagonable and denote a corresponding 
diagonaliser as M^, it follows that 

(2.30) Mo (I + pM(i))Mi •••(/ + p^+^i^f ') •••(/ + pKf^)Mk 

diagonalises A{p) for small p modulo 0{p^'^). 

Remark 2.5. If we had followed Remark l2.4l we could simplify this step (but with the cost 
of already determining the relevant matrices in the previous steps). If we always perform 
a perfect block-diagonalisation, the remainder term Aj^^^ has H^.i-block-structure and 
can be treated with one diagonalisation hierarchy instead of fc -|- 1 sub-steps. 
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2.2.3. A hierarchy of conditions and results. We recall that in each step of the iterative 

~ \k k] 

scheme there appeared one assumption, namely that A^+i = Af^^^ is diagonable. It is 
difficult to express these assumptions by conditions on the coefficients in the original 
expansion ()2.ip . We will abuse notation a little more and denote Aq = Aq. Then the 
following definition makes sense: 

Definition 2.2. The matrix family A{p) satisfies Assumption (A„) if the matrices A^ 
for k = 0,1, . . .n are all diagonable. 

Proposition 2.4. // A{p) satisfies (A„), then there exists a small interval In and a 
uniformly bounded family of invertible matrices Mn{p), p G In, with uniformly bounded 
inverse, such that 

(2.31) M-\p)A{p)Mn{p) 
is diagonal modulo 0{p^~^^). 

In Definition 12 . II we introduced the notion of a non-degenerate family of matrices. This 
definition can be relaxed by the following weaker notion. 

Definition 2.3. The matrix family A{p) is called non- degenerate of order n in /? = if 
it satisfies (A„), the matrix A„ consists of non-degenerate blocks and this number n is 
minimal. 

Obviously this definition coincides with the old one for n = 0. Furthermore, non- 
degeneracy of order n implies ( A„/ ) for all n' G N and allows for perfect diagonalisation 
by our scheme. 

Theorem 2.5. Assume that A{p) is non- degenerate of order n in p = 0. Then the 
statement of Theorem \2.1\ holds true. Furthermore, at least two eigenvalues of A{p) 
coincide modulo 0{p"') in p = 0. 

Remark 2.6. The statements of Propositions 12.21 and 12.31 transfer with the respective 
change in notation, so knowing an approximation of the diagonaliser M{p) allows us to 
estimate eigenvalues and eigenprojections of the matrix- valued function A[p). 

2.2.4. An example. We want to provide at least one detailed example to make the algo- 
rithm more comprehensible. For this we consider the matrix-valued function 

(2.32) A{p) = Aq + pAi + p^A2 
with coefficient matrices 

(2.33) Ao = -iaao\, ^i = -/? 

^ \0 0/ \6/2 5/2 

They will reappear in Section f3. 51 within the treatment of a model of thermo-elasticity. 

Step 0. The matrix Aq is degenerate with eigenvalues are 2a and twice. A diagonaliser 
of Aq is given by 

/l 1 0\ 1 /i 1 ^^ 

(2.34) Mo = 1 -1 , Mq-i = - 1 -1 

\0 1/ "^ \0 2, 



7\ 




/O 


7 , 


A2 = 





0/ 




\0 K 
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Hence 

(a 0\ /O /3 7 

(2.35) iM(p) =0 O+p/3 

\0 0/ \(5 

Step 1. For a 7^ we have Ho = (1,2) and the matrices Ai and K^ are given by 

(2.36) Ai = 0, Ki°l =1 

Transforming with / + pK\^ (and noting that Ai is aheady of diagonal form) yields 
iW(p) = (/ - rf + p\Kf^?)m{p){I + pKf) + 0{p') 
= Ao + p^(4nif]Ki°l) + 0(/,3) 

(a 0\ /^ 

(2.37) = +p2 -^ -^1 +0{p 



a 


U 


"\ 




H ^ I" 







u 








A 


+p' 





/32 


_ 


_/37 








0/ 


\ 





_^^ 


K 


7^ 



Step 2. Because Ai = we get Hi = Ho and the partition can not be refined. So the 
best we can do is to (1, 2)-block-diagonalise the matrix Ji[^J(/9) modulo 0{p^) in the first 
two sub-steps. Since it is already (l,2)-block-diagonal, we can skip these and proceed 

directly to the final one. The matrix A2 = ^2 ^ ^2 ^^ diagonable as soon as its lower 
right block is. So we calculate its eigenvalues. They are given by 

[2] _ P^ + jS 



(2.38a) X\ 



a 



and the latter two are distinct provided that (/3^ + {'j6 + ok))^ — Aan'yd / 0. We assume 
this; then assumption (A2) is satisfied and the matrix family A{p) is non-degenerate of 
order 2 in p = 0. 

It follows that eigenvalues and eigenprojections have full asymptotic expansions as p ^ 0, 
and the main terms can be read off from the above matrices, i.e., 

(2.39) Ai(p) = a + p^Afl + 0{p'), Xy,{p) = p^xfj^ + 0{p'). 

2.3. On the optimality of the conditions. We call a matrix family A(p) uniformly 
diagonable for all p if there exists a family of invertible matrices T(p), continuous in 
p, such that T~^{p)A{p)T{p) is diagonal and the matrices satisfy a uniform bound 

supp||r(p)|| < 00. 

Lemma 2.6. Assume that the matrix fam,ily A[p) is uniform,ly diagonable for all p < e 
up to and including p = 0. Then the assumptions (A„) are satisfied for all n. 
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Proof. If we plug in p = we obtain that Aq is diagonable and thus (Aq) follows. 
Assume now that for one particular n > the assumption (A„) holds, but (A„_|_i) fails. 
Then we can apply the first n iterations of our diagonalisation scheme and obtain the 
existence of a polynomial matrix function Mn{p), uniformly invertible for small p, such 
that 

n 

(2.40) M-\p)A{p)Mn{p) = Y, /^ A, + /,"+i4";ll + 0(p"+2) 

j=o 

is valid. Furthermore, ^^^i is not diagonable, i.e., one of its blocks is not diagonable. 
Note that on the subspace W corresponding to this block all diagonal matrices Aj are 
just multiples of the identity. 

By assumption the left hand side is uniformly diagonable and thus there exists a transfor- 
mation T{p) = M~^(p)T(p) diagonalising the left hand side. Restricting consideration 
to the above mentioned invariant subspace W gives the diagonal matrix 



(2.41) (Y^p^X^^)! + p-+'f-Hp)A^:fif{p) 



3 



+ (/>"+'). 
w 



Since the first addend is diagonal it remains to consider the last two. Dividing by /o"'^^ 
and taking the limit for p ^ (which exists due to continuity of T{p)) gives that ^^(O)!^ 
diagonalises the non-diagonable block and thus the desired contradiction. D 

Remark 2.7. From classical perturbation theory, see e.g. the book of Kato, [5], or Knopp, 
[6], it is clear that analyticity of A{p) implies that the eigenvalues are branches of alge- 
braic functions which have Puiseux series containing fractional exponents (with powers 
pk/p^ /c = 0, 1, . . ., where p corresponds to the size of irreducible groups of eigenvalues 
-so called A-groups- that are permuted if the degenerate point is encircled in the com- 
plex plane). If we assume that fractional powers do appear and no eigenvalues coincide 
identically (i.e., if we assume that the matrices are not permanently degenerate), then 
diagonalisers exist in a neighbourhood of p = 0, but due to the above lemma we know 
that they cannot be uniformly bounded / continuous in p. 

Remark 2.8. In Kato O Chapter II] asymptotic expansions of eigenvalues and eigen- 
projections corresponding to A-groups of eigenvalues are discussed and first terms are 
given. The approach used there differs from our treatment and is based on Dunford in- 
tegrals for resolvents and uses analytic dependence of the matrix family on the involved 
parameters. 



3. Applications 

Power stems from flexibility. We will show by a selection of applications to what extent 
the schemes introduced in Section [5] can be adapted to deal with real problems in the 
analysis of (partial) differential equations without loosing their constructiveness. The 
selection is not complete, but intended to give an impression of the variety of possible 
uses. 
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3.1. Hyperbolic polynomials. As a first application of the standard scheme of Sec- 
tion [2TT] we want to discuss the behaviour of the roots of hyperbolic polynomials for large 
spatial frequencies ^. 

Definition 3.1. A polynomial L{t, ^) in 1 + n variables and of degree m is called strictly 
hyperbolic if its m-homogeneous part Lm{T, ^) seen as polynomial in r parameterised by 
^ has 771 distinct real roots (j)i{^), . . . , (j)miO- 

We are interested in the behaviour of the roots of L{t,^), the so-called characteristic 
roots, as |^| -^ cxd. In particular, we want to prove the following asymptotic expansion: 

Theorem 3.1. The characteristic roots of a strictly hyperbolic polynomial have full 
asymptotic expansions as \S,\ — > oo, i.e., 

(3.1) r,{0 ~ mM + rfiv) + |^rV«(r?) + • • • , |^| - oo 

uniformly in rj = ^/|^| G §"~^ with smooth (algebraic) functions t- : S""-*^ — > C. 

Sketch of proof. The statement follows from the argument used to prove Theorem 12.11 
after rewriting it in terms of matrices and recognising that everything works uniformly 
in parameters. For this we rewrite the polynomial as 

m 

(3.2) i>(T,e) = c5]r— '=iciV(e), 

fc=0 

with \^\ Pk{0 polynomial in ^ of degree (at the most) k and po(C) = 1- Without loss of 
generalisation we set c = 1 and form a companion matrix with L[t, ^) as characteristic 
polynomial, 

/ 1 \ 



(3.3) £(e) = 1^1 



G 



1 

\-Pm{0 -Pm-l(0 -PliOj 

This matrix can be written as sum of homogeneous components 

(3.4) cio = lei A(r?) + £i(r?) + • • • + iei'^"/:™(^), V = m\ e s"~\ 

corresponding to the homogeneous parts of the polynomials \C\ Pk{C)- The assumption 
of strict hyperbolicity is equivalent to the fact that the matrix >Co(^) has m distinct real 
eigenvalues (/>! (r/) , . . . , (pmi'n) (which are uniformly separated by the compactness of S""^) 
and the standard diagonalisation scheme applied for p = |C|~^ gives representations for 
all functions involved in (13.11). D 



We see that the terms r ■ (rj) up to r ■ {rj) depend only on the homogeneous compo- 
nents Co up to Ck- Therefore, the following conclusion is apparent. 

Corollary 3.2. // the polynomial L{t,S,) — Lm{T,^) is of degree L < m — 1, then the 
coefficient functions r • (r/) up to r (rj) in the above given asymptotic expansion 

vanish identically. 
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The statements of Theorem 13.11 and Corollary 13.21 are of particular interest in the sit- 
uation of [Uj, where dispersive estimates for higher-order hyperbolic equations with 
constant coefficients are discussed. 

3.2. Asymptotic integration of systems of differential equations. Assume we 
are given a linear system of differential equations 

(3.5) v{t) = ^ = A{t)v, v{0) =voe C^, 

with a time-dependent coefficient matrix A{t) G C°°(M;C'"^™) having an asymptotic 
expansion 

(3.6) A{t) = AQ + t-^Ai + r'^A2 + ---, t^oo, 

with non-degenerate Aq. We assume further that derivatives of A{t) also have asymptotic 
expansions (which in consequence implies that we can differentiate the above expansion 
term by term). 

Theorem 3.3. Assume Aq is non- degenerate. There exists a uniformly bounded and 
invertible matrix function M{t) and a diagonal matrix function A(t), both having full 
asymptotic expansions as t ^ oo, such that the operator identity 

(3.7) (^_^(t))M(t) = M(t)(A_A(t)) mod f]0{t-'') 

holds modulo matrices decaying faster than all polynomials. 

Sketch of proof. We apply a variant of the standard diagonalisation scheme for diago- 
nalising this system. Thus, we construct recursively matrices 

(3.8a) Mkit) = Mo(/ + t'^Al'-^^ + t'^M^^) + • • • + t-^^M^), 

(3.8b) Afc(t) = Ao + t-^A(i) + i-2A(2) + . . . + t-^A^^), 

and follow the two steps from the proof of Theorem 12.11 Again Mq is the diagonaliser of 
Aq and Step 1 transfers directly. In Step 2 we include the differential operator in (j2.6p 

(3.9) B^>'\t) = (A _ A{t))Mk.i{t) - Mfc_i(t)(^ - Afc_i(t)) = 0{t-') 

and, thus, view it also as an operator identity. Analogousl}0, we define the matrices 
B(k) _ i[ia^^_^^t^MQ^B^^\t), A*^'^) = — diag-B^*-') and M*^'^) as solution to the commuta- 
tor equation [Ao,M''^''] = B^^' + A''^'. Since time-derivatives are one order better, the 
scheme works through, 

(3.10) - t-''{A{t)MoM^'''^ - MoM^^^Akit) - Mk^i{t)PS^^) 

- kt'^-^MoM^^^ 

is of order 0{t~^~^), and the recursion provides all matrices involved in the statement 
together with the necessary bounds. D 



Since we are concerned with ~A{t) in formula (|3.7|l we include a minus sign in the definition of A''^'. 



DIAGONALISATION SCHEMES AND APPLICATIONS 13 

Remark 3.1. If the matrix Aq is degenerate, multi-step schemes apply in a similar way. 

A particular application of Theorem 13.31 is the derivation of WKB approximations of 
solutions to hyperbolic systems. 



Corollary 3.4. Assume in addition that Aq is skew. Then the solutions to ()3.5p satisfy 

(3.11) v{t) = M{t) exp ( / A(s)ds J Q{t)M-^{t)vQ, 

with a uniformly hounded matrix Q{t) converging faster than polynomially to an invertible 
limit Qq = Q{oo) as t ^ cxd. 



Proof. We denote the remainder in ()3.7p as 

(3.12) R{t) = M-\t)M{t) - M^\t)A{t)M{t) + A(t) G p| 0{t-^) 



N 



such that v^^\t) = M-^{t)v{t) satisfies v'-'^\t) = {A{t) + R{t))v^'^\t). The diagonal part 
of this remaining system can be solved directly by means of the fundamental solution 

(3.13) £{t) = exp ( / A(s)ds J . 

To treat R{t) as perturbation we make an Ansatz for the fundamental solution of the 
diagonalised system of the form £{t)Q{t). This yields for Q{t) 

(3.14) ^Q(t) = n{t)Q{t) = {£-\t)R{t)£{t))Q{t), Q(0) = /, 
which can be solved in terms of the Peano-Baker series 

(3.15) Q(t) = / + V / 7^(tl) / 7^(t2)••• / " 7^(tfc)difc•••dt2dtl. 

i^^iJo Jo -/o 

If ^0 is skew it follows that £{t) and £^^{t) both satisfy polynomial bounds. Hence, 
TZ{t) decays fast in t, and the statement follows from the estimates 

(3.16) ||Q(t)||<exp('y 7^(s)dsV 

7^(s)ds exp i / 7^(s)dsj G f|0(t-^), 

in combination with Liouville's theorem, det(5(oo) = exp(J„ trace -R(s)ds) 7^ 0. D 

A variant of this approach was used in [231 Chapter 2] . For a slightly different diagonal- 
isation based method using less regularity of the coefficient see the treatise of Eastham, 
[H Chapter 1], or the second author's utilisation in |22j . 
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3.3. Diagonalisation within symbolic hierarchies. The method of the previous 
section can be extended to systems involving parameters. Typical applications are hy- 
perbolic partial differential equations with t-dependent coefficients, which are treated by 
means of the partial Fourier transform. In this case one has to be careful by choosing 
t large in dependence of the parameter and this leads to the introduction of so-called 
zones. We will not go into details here, but refer the reader to the fundamental treatise 
of Yagdjian on weakly hyperbolic problems, [231 Chapter 3], and applications deducing 
dispersive estimates for wave equations with variable propagation speed of Reissig and 
co-authors, [13], [9]. More involved considerations including several zones and different 
diagonalisation hierarchies turn up for the treatment of lower-order terms, e.g. in p], 
dU, [H] and [20] . 

A second modification of the standard scheme is to replace the matrices with complex 
entries by matrix-valued pseudo-differential operators and to apply it within certain 
symbolic hierarchies. This was developed in ^23j in order to prove well-posedness of 
hyperbolic systems with multiple characteristics, and similarly by Kubo- Reissig [8], [7], 
and Hirosawa-Reissig [3], to obtain corresponding results for strictly hyperbolic problems 
with certain non-Lipschitz coefficients. 

We will sketch the approach in a somewhat simplified case. We denote by Zhyp{N) and 
2pdiN) the subsets 

(3.18a) 2hypiN) = {(t,x,0 : i(0 > N}, 

(3.18b) Zpd{N) = {{t,x,0 ■■ m<N} 



of the extended phase space (0,T] x M" x M.?, where as usual (^) = y^l + |^p, and 
consider pseudo-differential operators corresponding to the following symbol classes: 

Definition 3.2. A symbol a{t,x,0 G C~((0,r] x M^ x R^) belongs to 5iv{mi,m2} if 
the symbolic estimates 

(3.19a) sup af5fafa(i,x,0 <Cfc,,,/3(e)'"^"l"l 7 

(3.19b) sup 9^"afa(t,x,0 <C;^(0"'' +""'"'"' 

it,x,^)eZpaiN) 

hold true for all multi-indices a, /3 S Ng and all k gNq. 

For a full account of calculus properties of such kinds of symbol classes we refer to j23| , 
we only mention embedding properties and relations to classical symbols here. Note that 
(0 > 1 and f-^ > T-^, and so 

(3.20a) Siy{mi — k, 1712 — i} ^^ Siyimi, 1712}, k,£>0. 

Furthermore, the definition of the hyperbolic zone implies 

(3.20b) S]\f{7ni — k, 1712 + k} ^^ Sn {1^1,1112}, k>0 

The embedding hierarchy (I3.20ap with i = will be denoted as calculus hierarchy. It 
allows us to transfer the usual symbolic calculus and corresponding asymptotic expan- 
sions from the Hormander classes 5™q(M^ x M?) via the embedding 5iv{mi,m2} ^-> 
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C°°((0, T], 5™Q (R" X M?)). Thus, the composition of operators corresponds to the usual 
Leibniz product (t of symbols. Symbols from S]sf{mi,m2} which are invertible modulo 
smoothing operators from S]\f{— 00,1712} with a parametrix in S]s[{—mi, —7712} will be 
called elliptic, and the parametrix is denoted with an upper ft- The second hierarchy will 
be used within the diagonalisation scheme, and we denote the corresponding residual 
class of this hierarchy by TIn{^ = flfc '^Arjm — k,k}. 

Assume for the following that we have a matrix- valued symbol A{t,x,S^) £ Sn{^,0}, 
whose eigenvalues satisfy the non-degeneracy assumption 

(3.21) \Xi{t,x,0-Xjit,x,0\>C{0 
uniformly in i ^ j and {t,x,(,) € Zhyp{N) for a certain large A^. 

Theorem 3.5. There exists an elliptic matrix valued symbol M(t, x, ^) G iSatJO, 0} and 
a diagonal symbol A(t, x,^) G 5Ar{l,0} such that the operator identity 

(3.22) (A_^(t,x,e))HM(t,x,0=M(t,x,OS(^-A(t,x,e)) mod W^{1} 
holds modulo symbols from the residual class WatJI}. 

Sketch of proof . The proof transfers almost word by word from Theorems l2.1l and l3.3l ex- 
cept that we now replace multiplications by the Leibniz product [j and solve the Sylvester 
equation ()2.7p only modulo symbols of lower order (such that (12. 8p still remains true) . 

For completeness we give the main steps and the corresponding symbol estimates. We 
construct 

(3.23a) Mk{t,x,0 = Mo{t,x,odl + Y.M^J\t,x,0 

and 

k 

(3.23b) Afc(t,x,0=Ao(t,x,0 + J^A(^)(t,x,e) 

i=i 

with matrix-valued symbols satisfying Mq G 5Ar{0,0} elliptic, Aq G 5Ar{l,0}, M^^' G 
Sn{-JJ} and A(^) G Sn{1 - j,j}. 

In Step 1 we choose MQ{t,x,^) to be a diagonaliser of the full symbol A{t,x,^) within 
Zhyp{N), which is uniformly bounded and invertible. This can be done by (|3.2ip . Fur- 
thermore, we define Aq = diagMgd^tt-^o- In Step 2 we proceed recursively and, in 
analogy to ()2.6p . define for A; = 1, 2, . . . 

(3.24) B'^''\t,x,0 = {^^-A{t,x,0)mk-i{t,x,0-Mk-i{t,x,m{^^-^k-i{t,x,0) 

For fc = 1 we know from Step 1 that B^^' G 5{0, 1}. For our recursive argument we 
assume S^'^) = M^^^B^^^ G Sn{1 - k,k} and set A('=) = -diag^('=). Defining M^^^ by 
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(12:81) 



{(B('=){t,x,0) 
0, i = j, 

with XN(t,x,(,) a smooth cut-off function locahsing to Zhyp{N), imphes 

[Ao, M^'^)] = AoHM^'^) - M^'^^SAo 

(3.26) =^W+A('=) mod5jv{-A:,/fc}, 
so 

^(fc+i) ^ Qik) _ A^Mo^M^''^ + Mo'^M^'^HAo + MoHA(^) 

(3.27) e Mott (^('^^ - [Ao, M^'^)] + A^'^)) + SN{-k, k + l} 

eSN{-k,k + l}. 

We can proceed inductively to construct all M^, A^ and the statement itself follows by 
forming asymptotic sums instead of the finite ones in (j3.23p . D 

This statement may be used to deduce well-posedness of degenerate Cauchy prob- 
lems. Assume that the matrix A{t,x,^) £ Sm{^,0} is hyperbolic in the sense that 
\ReXj{t,x,^)\ < C uniform in (t,x,^) G (0,r] x R^ x M^ in combination with (fOT]) . 
but degenerates as t approaches 0, 

(3.28) \\A{t,x,0\\<l, \\dtA{t,x,m^]- 
We may ask whether we can pose the Cauchy problem 

(3.29) ^v(t,x)=A(t,x,B)v(t,x), w(0, •) G L^(M"; C"). 
dt 

The standard energy argument does not work, because this would mean that we have to 
differentiate the full symbol A{t,x,^). But the above diagonalisation argument simpli- 
fies the problems substantially. If we consider v^^'{t,x) = M\t,x,D)v{t,x) and apply 
Theorem 13. 5t we obtain 

(3.30) —v^^\t,x) = {A{t,x,B)+R{t,x,B))v^^\t,x), 

with a remainder R{t,x,S,) G TCn{^}- The diagonal part of this system can be solved 
in terms of a diagonal matrix of (elliptic) Fourier integral operators of order zero. This 
follows from the hyperbolicity assumption. Furthermore, it can be shown that TC]s[{m} is 
invariant under conjugations with such Fourier integral operators (see |23t Prop. 3.8.12]), 
and for an arbitrary P G TIn{^} the operator ^—P{t, x, D) has a pseudo-differential fun- 
damental solution in C([0, T], ^°(M")), included (see ^ Prop. 3.9.1]). Well-posedness 
follows. 

The argument can be found, in full detail and including fast oscillations within the 
coefficient (i.e., a further log-term in the estimates), in [8]. 
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3.4. Pseudo-differential decoupling of hyperbolic-parabolic coupled systems. 

The two-step diagonalisation scheme appeared first in Wang |17| . |18j in appHcations 
decoupling hyperbolic-parabolic coupled systems, and was used to study the propagation 
of singularities, [10]. In [12], [21] it was an essential tool for understanding dispersive 
estimates for anistropic thermo-elasticity in two space dimensions. 

Again we will give a somewhat simplified example, which shows the main arguments 
involved and the importance of the two-step procedure. For this we choose a thermo- 
elastic system, 

(3.31a) Utt - {2fi + A)VV • [/ + /iV x (V x [/) + 71 V6' = 0, 

(3.31b) Ot - (i^M + 72V • C/t = 0, 

with smooth scalar functions A, /i,/3, 71,72 G C°^(Ri x Q,), Q, c M^, subject to the 
restrictions // > 0, 2/i + A > 0, /3 > and 7172 > 0. If one employs a Helmholtz 
decomposition of the displacement vector U , 

(3.32) U = UP + U', V xUP = 0, V • [/' = 0, 

either globally on R'^ or locally on any small open subset Q' of the underlying domain 
ri with trivial second cohomology class, the system simplifies. The solenoidal part U^ 
satisfies the (well-studied) wave equation C//j — fiAU^ = 0, while the potential part 
UP = VV satisfies 

(3.33a) VVtt - (2/i + X)VAV + 71 V6' = 0, 

(3.33b) Ot - /3^Ae + -f2AVt = 0. 

We can write this as a first order 3x3 pseudo-differential system of a particular structure. 
We choose for this as new unknown W = {\/—AVt, AV, O)"^ such that 

TZ~\2fi + X)V 7^~l7lV^ 
(3.34) ^W = \ -^/^A \W = A{t,x, !})]¥, 

p^A 




with n : ^/^AV ^ VV the/a Riesz transform (cf. [13 Chapter III] for n = M" or 
use any invertible pseudo-differential operator with principal symbol S,/\(,\ if O 7^ M"'). 
The symbol A(t, x, ^) belongs to the class C°°(]R, Siq(T*Q')), but the only second order 
entry is at the lower right corner. 

An adaptation of the two-step scheme of [T7] allows one to construct an elliptic matrix- 
valued pseudo-differential operator M G C°°(M, S^qIT*^'))) such that 

(3.35) r^ - Ait, X, o) m{t, X, = M{t, X, m (^ - Mt, X, 

holds modulo C°°{R, S'{;^(T*f^'))> where Ai = b-diag(^ii, 012) consists of a 2 x 2-block 
All G C~(M,5{o(T*0'j) and a scalar entry 012 G C°°(M, 5f o(T*f7'))- The 2 x 2 block 
is strictly hyperbolic, and can be studied further by diagonalisation, i.e., there exists a 
second diagonaliser N{t,x,(,) G C~(M, 5?o(T*fi')) such that 

(3.36) ("^ - Ai{t, X, 0) m{t, X, = Nit, x, e)« f^ - A2it, x, 
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holds modulo C°°(]R, S'j~^(T*il')), with a diagonal matrix A2 = diag(a2i, 022, 012) having 
entries 0,21 G C°°(M, S'|Q(T*il')), i = 1,2, and ai2 as above. Knowing the principle part 
of these scalar symbols allows to describe the propagation of singularities. 

3.5. Dispersive estimates and diffusive structure for thermo-elastic models. 

Multi-step schemes have been developed in ^j in order to study models of thermo- 
elasticity with additional terms of lower order. They have been used in [12], [21] and 
will also be central to the understanding of the generalisation of these results to higher 
dimensions. 

Again we will provide an application. It is taken from [5] and concerned with the 
derivation of dispersive estimates for one-dimensional thermo-elastic systems with terms 
of lower order. We restrict ourselves to one model, classical one-dimensional thermo- 
elasticity with an additional damping term: 

(3.37a) utt - t'^Uxx + liOx + mut = 0, 

(3.37b) 6t - K.9r,x + 72Utx = 0, 

together with initial conditions u(0, •) = ui, ut{0,-) = U2 and ^(0,-) = ^o- Parameters 
T, K, 7i, 72 and m are assumed to be positive constants. Applying a partial Fourier 
transform with respect to the x-variable reduces the system to a system of ordinary 
differential equations, which can be written as system of first order mV = {u-^-,u^, 9)'^ , 
u± = ut i: iT^u. A short calculation gives 

(3.38) ^V = A{S)V={A^ + iAi+eA2)V, 



with matrices 




(3.39) 






-m 0' 


-m 




















— K 



^1 = i I -T 71 I , ^9 

72 

If we assume for simplicity that A{^) has no multiple eigenvalues, we can represent 
solutions to this system as a sum 

(3.40) V{t,i)= Y. ^'"^^^PhO^o 

!^(?)especA(5) 

over the spectrum of the matrix A{^). In order to understand properties of solutions such 
as dispersive estimates or descriptions of asymptotic profiles it is enough to calculate 
the eigenvalues z^(C), or at least to describe their main properties. 

This can be done in three steps. First, a consideration of the characteristic polynomial 
of A{S^) implies that no purely imaginary eigenvalues occur for ^ 7^ 0. Thus solutions 
with bounded frequencies away from either decay exponentially or increase exponen- 
tially. In a second step one can diagonalise and determine asymptotic expansions of 
eigenvalues (and eigenprojections) as ^ — > and ^ — > 00. This proves that Rei^(^) < 
everywhere and gives enough information about the behaviour of the eigenvalues as 
^ ^ to determine decay estimates in a third step. 
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We will not give the precise calculations, for which see [4J. It is easily seen that A{^) is 
non-degenerate of order 1 as ^ — > oo, and, following the first two steps of the procedure 
from Section 12. 2^ the eigenvalues have asymptotic expansions of the form 

(3.41) ,^par{0 = -^e + ^ + O(r'), 

(3.42) i^hypAi) = ±i^^ - ^ + 0(r'). 

Thus, their real part is uniformly negative for large ^, which implies exponential decay 
for large frequencies. It remains to consider ^ — > 0. This calculation has already been 
done in Section 12. 2. 4^ where it was also seen that A(^) is non-degenerate of order 2 
(under the above restrictions on the parameters). It follows that the three eigenvalues 
satisfy 

(3.43) y^[^) = -m + \oe + 0{e), y±{i) = -\±f + 0{e), 

with positive constants Aq, X± depending on the given parameters. Thus, one mode 
also leads to exponential decay, while the other two resemble a parabolic-type behaviour 
close to the corresponding heat equations vt = X±Vxx- 

The above mentioned properties of the eigenvalues of A{p) in combination with bound- 
edness properties of the Fourier transform and Holder inequality imply directly the 
following a-priori estimate. Again the detailed proof can be found in [3]. 

Theorem 3.6. The solutions to the the system (j3.37p satisfy the a-priori estimate 

(3.44) \\{uuUx,e){tr)\U<C{l + t)-'^'^'p'-^\\\ul\\Hv,r+^ + \\u2,e4Hv,r) 

for all indices l<p<2<q<oo and r > {1/p — ^/q)- 
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